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Screening properties of the quark gluon plasma are studied from the Polyakov-Hne correlations in 
lattice QCD simulations with two flavors of improved Wilson quarks at temperatures T/Tpc ~ 1-4 
where Tpc is the pseudocritical temperature. Using the Euclidean-time reflection symmetry and 
the charge conjugation symmetry, we introduce various types of Polyakov-line correlations and 
extract screening masses in magnetic and electric sectors. We find that a ratio of the screening 
masses in the electric sector to the magnetic sector shows qualitative agreement with a prediction 
from the dimensionally-reduced effective field theory and the ,y\^ — 4 supersymmetric Yang-Mills 
theoryatl.3<r/rpc<3. 
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1. Introduction 

Electric and magnetic screening masses of the gluon are the fundamental quantities of the 
quark-gluon plasma at high temperature In this report, we employ lattice QCD simulations with 
two flavors of improved Wilson quarks at temperatures T/Tpc — 1-4 and try to extract both screen- 
ing masses from the correlations of the Polyakov-line operator classified under the the Euclidean- 
time reflection and the charge conjugation Also, we make a comparison of our results with the 
predictions from the dimensionally-reduced effective field theory and the ^ = 4 supersymmetric 
Yang-Mills theory. 

2. Correlations of the Polyakov-line operator 

In order to extract the electric and magnetic screening masses, we classify the Polyakov-line 
operator into different classes on the basis of the symmetries under the Euclidean time-reflection 

and the charge conjugation (C). Both are good symmetries of QCD at zero chemical potential. 
(Note that corresponds to the product of the time reversal T and charge conjugation C in the 
Minkowski space.) Under J^, the gluon fields are transformed as, 

A(t,x)^A(-t,x), A4(t,x)^-A4(-t,x). (2.1) 

An operator belongs to the magnetic (electric) sector if it is even (odd) under On the other 
hand, the charge conjugation changes the gluon field as, 

A^(t,x)^-A;(t,x). (2.2) 

The Polyakov-line operator 

"l/T 



n(x) = Pexp 



ig / JtA4(t,x) 

JO 



(2.3) 



transforms under ^ symmetry and C symmetry as 

^ : D.^D.\ C : D.^D.*. (2.4) 
Therefore, we define =^-even (odd) operator by its Hermitian (anti-Hermitian) component: 

aM = ^(n + al'), aE = ^(a-a"0, (2-5) 

where the subscript implies that the operator is in the magnetic (electric) sector. Furthermore, they 
can be decomposed into C-even and C-odd operators as 

aM.± = \{^M±^*M), I1e,± = ^(^^e±^^e)- (2-6) 



We can define two types of gauge invariant correlators using Eq.(2.6) as 



CM,+ (^,r) = (TrnM.+ (x)TrQM,+(y)) - (Trn)^, (2.7) 
4_(r,r) = (TraE,-(x)TraE,-(y)), (2.8) 
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where r = |x — y|, and the superscript "I" impUes that the correlators are gauge invariant. Since 
Tri^M - = TrnE.+ = 0, the correlators of these operators are identically zero. If we are allowed to 
fix gauge, we can construct four types of correlators: 

CmA^'^T) = (TrQM,+ (x)nM,+ (y)) - (Tra)2, (2.9) 

CM,-ir,T) = (TrnM,-(x)nM.-(y)), (2.10) 

Cl4r,T) = (TraE,+(x)nE,+ (y)), (2.11) 

Cl4r,T) = (TrnE,-(x)aE,-(y)). (2.12) 

where the superscript "F" implies gauge-fixed correlators. In this study, we adopt the Coulomb 
gauge fixing. 

We define the screening masses of these correlators through the fit with the following screened 
Coulomb form at long distances: 

ci%T) ^ ajf) (D^— (2.13) 

where a^^^(r) and m^^(r) are fitting parameters. Then we can extract six screening masses all 
together, 

(Magnetic sector) : m\^ _^_, m^^, (2.14) 
m'^^ (Electric sector) : m^ _, m| _|_, m^ _. (2.15) 

3. Numerical simulations 

We employ a renormalization group improved gauge action and a clover improved Wilson 
quark action with two flavors. The simulations are performed on a lattice with a size of A^,^ xA^, = 
16-' X 4 along lines of constant physics, i.e. lines of constant mps/my (the ratio of pseudoscalar 
and vector meson masses) at T = in the space of simulation parameters. Details of the lines of 
constant physics with the same actions are summarized in Refs. ^, ^. We take two values of 
mps/mv = 0.65 and 0.80 with the temperature range of T /Tpc 1.0-4.0 (9 points) and 1.0-3.0 (7 
points), respectively, where Tpc is the pseudocritical temperature along the line of constant physics. 
The number of trajectories for each run after thermalization is 5000-6000, and we measure physical 
quantities at every 10 trajectories. The fits are performed by minimizing X^/^df with fit ranges of 
0.35 < rr < 1.0 for _ and 0.5 <rT < 1.25 for other orrelators. 

Left (right) panel of Fig. |l] shows results of screening masses in magnetic (electric) sector at 
mps//Mv = 0.65 as a function of temperature. As temperature increases, the magnetic screening 
masses converge to a similar value, whereas the electric screening masses shows wig ^ ~ nig _ < 
wig . In the weak coupling limit, these properties may be interpreted as follows. The Polyakov-line 
operator in the magnetic (electric) sector has even (odd) powers of g: 

nu l + ^(^igj^^^^ dTA4^ +■■■, (3.1) 
-ig j^J^ dxA^-jlig j^J^ dxA^ +•••. (3.2) 
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Figure 1: Screening masses in magnetic (left) 
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electric (right) sectors at mps/my = 0.65. 



Then all the Polyakov-line correlations in the magnetic sector start with the correlation of the 
second powers of A4: 

CUrJ) - /(TrA2(x)A2(y)) + 0(/), (3.3) 
Cl{r,T) oc g\TrAl{x)TrAl{y)) + 0{g'). (3.4) 

On the other hand, Polyakov-line correlations in the electric sector have different powers depending 
on the gauge choice, because T1A4 = 0: 

Cl{r,T) oc ^2^TrA4(x)A4(y)) + (3.5) 
Ci{r,T) oc /(TrAl(x)TrA3(y)) + (3.6) 

Then the leading contribution in Cg comes from single temporal-gluon exchange, whereas that in 
Cg comes from triple temporal-gluon exchange. This may lead to the exponential fall of Cg _ much 
faster than Cg ^ and hence the mass ordering in Fig.[T| Further discussions on a relation between the 
Polyakov-line correlations and the gluon propagators in the weak coupling are given in Appendix. 

3.1 Comparison with 3D effective field theory and ,jV = A supersymmetric Yang-Mills 
theory 

Let us compare our results with the predictions by the dimensionally reduced effective field 
theory (3D-EFT) and the ^ = 4 supersymmetric Yang-Mills theory (SYM). Here we focus our 
attention on the screening masses with (=v^,C) = (+,+) and (—,—), i.e. myi,+ and mg.-, which 
are channels calculable from the gauge invariant Polyakov-line correlations ^ In the 3D-EFT ap- 
proach the screening masses at T = 2A^ with being the QCD scale parameter have been 
calculated for Nf = 2 as 

3D-EFT {Nf = 2) : mM,+ ~ 3.60r, me - ~ 6.467. (3.7) 

'Although our Polyakov-hne operators do not have definite angular momentum J and the parity P, we assume 
that the asymptotic behavior of our Polyakov-line correlations picks up the ground state contribution corresponding to 

(y,p) = (o,+) 
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Figure 2: Comparison of the screening ratio, '«e,-/wm,+, with predictions in the dimensionally-reduced 
effective field theory (3D-EFT) [|] and ^ = 4 supersymmeti-ic Yang-Mills theory (SYM) 



The screening masses of the Polyakov-line correlations in = 4 supersymmetric Yang-Mills 
theory in the limit of large Nc and large 'tHooft coupling (A = g^N^) were calculated by AdS/CFT 
correspondence [^: 

^ = 4 SYM : mM,+ ~ 7.34r, wie,- ~ 16.ir. (3.8) 

Our results of the screening masses shown in Fig. [l] indicates that 

4D-lattice QCD {Nf = 2) : mM,+ =AT- 6T, me.- = lOT - 14T. (3.9) 

In all three cases, we have the inequality: mM,+ < m^,-- Since we cannot compare the absolute 
magnitude of the screening masses in QCD and those in .J/' = 4 SYM because of the different 
number of degrees of freedom, we take the ratio, mE-/mM.+ in Fig. § for above three cases. We 
find that, in spite of the difference in the magnitude, the screening ratios agree well with each other 

for 1.5 <r/rpc <3. 
4. Summary 

We investigated the screening masses in hot quark-gluon plasma from the Polyakov-line cor- 
relations classified by the Euclidean-time reflection (J7) symmetry and the charge conjugation (C) 
symmetry. The screening masses are calculated in the magnetic (J^-even) and electric (^-odd) 
sectors from the lattice simulations of Nf = 2 QCD above Tpc with the RG-improved gluon action 
and the clover-improved Wilson quark action on a 16^ x 4 lattice. We found that, the magnetic 
screening masses obtained by the gauge-invariant and gauge-fixed correlations coincide with each 
other for T /Tpc > 2, whereas the electric screening mass obtained from the gauge invariant cor- 
relator is larger than that obtained from the gauge fixed correlator. Our results obtained from the 
gauge invariant correlators are compared with the calculations in the dimensionally-reduced effec- 
tive field theory [^] and the ^ = 4 supersymmetric Yang-Mills theory [^]. We found that the ratio 
of the electric and magnetic screening masses in three cases are qualitatively consistent with each 
other for at 1.3 < T /Tpc < 3. Further discussions between the Polyakov-line correlations and the 
gluon propagators are given in Appendix. 
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A. Appendix: Relation to gluon propagators 

In this Appendix, we consider a direct relation between the Polyakov-line correlations and the 
screening mass in gluon propagators. Here we define the magnetic (electric) screening mass, wjm 
(niE), from the spatial (temporal) gluon propagator at asymptotic spatial distance, 

{Mx)My)) ^ T (|x-y| =r^oo), (A.l) 

(A4(x)A4(y)) ^ r ^ (|x-y|=r^oo). (A.2) 

In the weak coupling limit, the gauge fixed correlation C| _ is related to the electric propagator as 

Cl_{rJ)^a{T)—^. (A.3) 

On the other hand, the gauge invariant correlation Cj^ _|_ would have not only an exchange of two 
temporal gluons but also an exchange of two spatial gluons : 

Although b2 has higher powers in g, the exchange of the spatial gluons could dominate at long 
distances if mm < rriE- 

We extract by fitting Cg _ to the right-hand-side of Eq. (A3) with the fitting parameters, 
a{T) and mE{T). Then we extract tum by fitting a combination of the Cg _ and as, 

:ci+C2(r)exp[2(m£(r)-mM(r))r], (A.5) 



(cL(r,r))2 



with ci, C2 and rriM being the fitting parameters. The fit ranges are chosen to be 0.5 < < 1.25. 

Results of m£(r) and mM{T) are shown in Fig. || (left) as a function of T /T^ at mps/my = 
0.65. Similar results are obtained for wips / my = 0. 80 too. We find that mm is smaller than at all 
temperatures we calculate, which is consistent with a prediction in the thermal perturbation theory: 
Mm ~ 0{g^T) < rriE ~ 0{gT). We also find that, for T /T^^ > 2, both screening masses decreases 
as T increases, whereas me and mm behaves differently for 1 < T/Jpc < 2, i.e. me increases and 
niM decreases when the temperature approaches to Tpc- The behavior of mm is not expected from 
the leading-order perturbation theory. 

In the right panel of Fig. ^, me and niM calculated directly from the gluon propagators in the 
quenched approximation are recapitulated Although niM < niE is also found for T /Tc ^ 1.2, 
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Figure 3: Results of the electric and magnetic screening masses calculated from the Polyakov-line correla- 
tions in Nf — 2 lattice QCD at mps/mv =0.65 (left) and from the gluon propagators in Nf — lattice QCD 
[i (right). 



tnM ('«£:)increases (decreases) as T approaches to in the quenched calculation. Such qualitative 
difference between Nf = 2 and Nf = near the (pseudo-) critical temperature may be attributed to 
the different behavior of the the phase transition in two cases. 
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